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We find an analytical expression for the conductance of a single electron transistor in the regime 
when temperature, level spacing, and charging energy of a grain are all of the same order. We 
consider the model of equidistant energy levels in a grain in the sequential tunneling approximation. 
In the case of spinless electrons our theory describes transport through a dot in the quantum Hall 
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the peak's maximum as a function of temperature, and the values of the conductance in the odd 
and even valleys. 
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I. INTRODUCTION 

Recent progress in mesoscopic fabrication techniques 

I has made possible not only the creation of more sophis- 
ticated devices but also greater control over their prop- 

I erties. Electron systems confined to small space regions, 
quantum dots, and especially their transport properties 
have been studied extensively for the last decade;-'^ In 
particular, an individual ultra-small metallic grain of ra- 
dius less than 5 nm was attached to two leads via ox- 
ide tunnel barriers, thus forming a single electron tran- 
sistor (SET);'^ Applying bias voltage, V, between two 

I leads allows one to study transport properties of the sys- 
tem. Fig. n Alternatively, a SET can be formed by de- 
pleting two-dimensional electron gas at the interface of 
GaAs/AlGaAs heterostructure by applying negative volt- 
ages to the metallic surface gates. ^ 

In this paper we will assume that the bias voltage is 
infinitesimally small, V ^ 0. This corresponds to the 
linear response regime. In order to tunnel onto the quan- 
tum dot, an electron in the left lead has to overcome a 

\ charging energy, Ec = e^/2C, where e > is the el- 
ementary charge; C is the capacitance of the quantum 

' dot. If T <C Ec then conductance through the system is 
exponentially suppressed. This phenomenon is called the 




FIG. 1: Scheme of the Coulomb blockade setup. 



Coulomb blockade. However if we apply a voltage, Vg, to 
the additional gate capacitively coupled to the dot, the 
Coulomb blockade can be lifted. Indeed, changing Vg one 
can shift the position of energy minimum so that ener- 
gies of the quantum dot with and N^ + l electrons will 
become equal and an electron can freely jump from the 
left lead onto the dot and then jump out into the other 
lead. Thus, current event has occurred and a peak in 
the conductance, G, corresponding to this gate voltage is 
observed. By changing the gate voltage one can observe 
an oscillation of the conductance or Coulomb blockade 
oscillations. 

One-particle energy levels in the quantum dot, {Ei}, 
are given by the solution of the Schrodinger equation in 
the quantum dot's potential. The mean spacing between 
these energy levels is SE. The conventional assumption 
that Ec ^ SE is not valid in the case of sufficiently small 
dots. In fact, in the recent experiments?^ where a Cgo 
molecule has acted as a quantum dot, the level spacing 
is of order charging energy. Experimenlji was performed 
at T = A.2K as well as at room temperature. In other 
experiments^'^ with quantum dot formed by depleting 
2DEG^ and ropes of carbon nanotubes acting as a quan- 
tum dot,^ charging energy is only three times larger than 
the spacing SE. 

Though Coulomb blockade oscillations have been stud- 
ied in a number of important limiting cases ^iSs^iifl the 
problem in the case when values of Ec, SE, and T are 
all of the same order has not been theoretically addressed. 
Let us note that energy levels of the quantum dots are 
random and obey Wigner-Dyson statistics with the fluc- 
tuation of order of their mean.^^ Nonetheless to go as far 
as possible in the analytical treatment of the problem we 
have to assume that energy levels in the quantum dot 
are equidistant. In this paper we derive an analytical 
expression for the linear conductance, G — I/V\v^o, in 
the case of spinless as well as spin-i fermions. 

In Sec. II, we describe our model and the assumptions 
involved. We write the model assuming spin-i fermions. 
In Sec. Ill, we consider the linear conductance in the case 
of spinless fermions. We obtain an analytical expression 
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for the conductance and analyze its limiting cases. In 
Sec. IV, we consider one possible application of the Sec. 
Ill results, namely tunneling through the edge states in 
a quantum dot placed into a strong magnetic field. In 
Sec. V, the linear conductance as well as its properties 
in the case of spin-i fermions is considered. In Sec. VI, 
we summarize our findings. 



II. THE MODEL 

Hamiltonian of the system in question is 

H ^ Hi+Hd + f. (1) 

Here, the first term is the Hamiltonian of noninteracting 
electrons in the left and right leads: 

Hi^Yl ^''^-^Cfc^ + Epcl„Cp„, (2) 

kcr per 

where a continuum of states in each lead, \ka), \pa) is as- 
sumed; Ek, Ep and Cka, Cpa are the energies and electron 
annihilation operators in the left and right leads, respec- 
tively; a stands for the ^-component of spin. The chem- 
ical potentials of the leads, /i ^ Ec,SE,T, are shifted 
according to the bias voltage, V, applied, Fig. [3 We will 
assume that leads are in thermal equilibrium at tempera- 
ture T and, thus, occupied according to the Fermi-Dirac 
distribution. 

The second term in Q is the Hamiltonian of the quan- 
tum dot: 

Hd^J2 ^'4c»- + U, (3) 

where first term is the kinetic energy of electrons in the 
quantum dot: {Ei} is a discrete set of the quantum 
dot's energy levels; Ci^'s are the annihilation operators. 
The second term, U describes the electron-electron in- 
teraction in the quantum dot. We adopt the simplest 
model for the interaction, namely, the constant interac- 
tion model. In this model the Coulomb interaction of the 



electrons depends only on the total number of electrons 
in the quantum dot: 

U{N) = EcN^ - eVeN, (4) 

where N = cL'-if ~ total number of excess 

electrons; Ni is the total number of positively charged 
ions. The second term is the contribution from external 
charges. They are supplied by the ionized donors and 
the gate: Ve = Vd + aVg, where a is a function of the 
capacitance matrix elements of the system. Thus, Ve 
can be varied continuously by changing gate voltage, Vg. 
U{N) can be rewritten as 

U{N) ^ EciN - Ngf + Const, (5) 

where Ng = eVe/2Ec is the dimensionless gate voltage. 
The third term in is the tunneling Hamiltonian: 

T = ^ {tkic\„Cic, + h.c}j + ^ {tpicl^Cia + h.c.) , (6) 

where tki and tpi arc matrix elements of tunneling into 
the left and right leads, respectively. 

We assume that the dot is weakly coupled to the leads; 
that is, the conductances of the dot-lead junctions are 
small: G'"'' <C e'^ /h, where h is Planck's constant. Equiv- 
alently, the widths of the quantum dot's energy levels 
contributing to the conductance, F^ = F' -I- F[, must 
be small compared to spacing between them: F^ <C SE. 
This, together with Fj ^ T assumption, allows us to 
characterize the state of the dot by a set of occupation 
numbers, {riia-}^ 

III. LINEAR CONDUCTANCE IN THE 
SPINLESS CASE 

The model formulated above has been studied by 
Beenakker in the sequential tunneling approximation^ifi 
that is, conservation of energy was assumed in each tun- 
neling process, and cotunneling was neglected. There- 
fore, to find the stationary current, kinetic equation con- 
siderations can be applied. In the linear response regime 
an analytical formula for the conductance has been ob- 
tained. In the case of spinless fermions>ifi 
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x[l-nF{E,-fi + U{N)-U{N~l)], (7) 

where 

F^ - 2^ ^ \tk,\^5 [E, -Ek + U{N) - U{N - 1)] 

k 

and 



FIG. 2: Electrostatic potential energy along a line through 
the tunnel junctions. 



F[ = 2TrY,MS [E, ~Ep + U{N) ~ U{N ~ 1)] 
p 



3 



are widths of the quantum dot's level i associated with 
tunneling into the left and right leads, respectively; 
Peq{Ne) IS the equilibrium probability that the quantum 
dot contains electrons; Ff,q{Ei\Nf,) is the occupation 
number of level i given that the dot contains iVg elec- 
trons; npiE) is the Fermi-Dirac distribution; and /i is 
the chemical potential in the leads. 

The quantity Feq{Ei\Ne) in O is the most non-trivial 
one to calculate. It is the occupation number of the level 
i in the canonical ensemble (iVe is fixed). In the limit 
SE/T 0, Feq{Ei\Ne) bccomes a Fermi-Dirac distribu- 
tion with the appropriately chosen chemical potential: 
p, = {Eq + Ei)/2, where Eq corresponds to the energy of 
the last occupied energy level at T = 0, Fig.Ot; Ei cor- 
responds to the energy of the first empty energy level at 
T = 0. In the opposite limit 5E/T oo, the Fermi-Dirac 
distribution with fl = {Eq + Ei)/2 apparently breaks 
down: the occupation number of level j = 1, for exam- 
ple, see Fig. Et, is m = e'^^/'^, not e'''^/^^ as the 
Fermi-Dirac distribution would predict^ 

The occupation number in question isi^ 
FeqiE,\N,) = ^_^P,,({n,})5„.,i5^^^^„^ 

{"J 



where Ppq[{ni\) is the equilibrium probability of the 
state of the quantum dot; (3 = l/T; and the de- 
tailed definition of F{Ne) will follow. The reason for 
writing this equation is to show that analytical calcu- 
lation of the occupation numbers is hardly possible for 
arbitrary quantum dot's energy level structure, {Ei}. 

The only way to overcome this difficulty is to assume 
that energy levels in the quantum dot are equidistant. 
Then one can use the bosonization techniqua^^ (see Ap- 
pendix) to find the exact analytical expression for the 
occupation numbers in the canonical ensemble. It was 



done by Denton, Muhlschlegel, and Scalapino 
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FIG. 3: a) occupation numbers in the canonical ensemble for 
the equidistant energy levels at T = 0; b) mapping sum over 
i onto sum over j {N = 2 case is shown) . 



r[m^ + (2j-l)ml4f 



where 



Ei-C 1 

__ + _ = integer. 



(9) 



(10) 



5E/T = S = const, C is the energy corresponding to the 
highest occupied energy level at T = plus SE/2, Fig.|3i. 
This quantity C is somewhat similar to the chemical po- 
tential of a dot, though, strictly speaking, the chemical 
potential is not well-defined for a dot in the canonical 
ensemble. The difference, C ~ is a linear function of 
the gate voltage. Therefore, by properly adjusting "zero" 
value of the gate voltage it can put to zero. Hereinafter, 
we assume that C — = 0. 

To calculate the conductance we also need to know 
PeqiNe), the probability that the dot, in thermodynamic 
equilibrium with the reservoirs, contains electrons. It 
can be calculated in the grand canonical ensemble: 



Peq{N,) = -^e 



(11) 



where = J2n =o e"'^'^^^"^ is the grand partition func- 
tion; (p{Ne) is the thermodynamic potential of the quan- 
tum dot. It can be expressed via free energy of the dot's 
internal degrees of freedom, F{Ne): 



hence, 



where 



Z{N, 



Peq{N,) = —e 



-I3U{N) 



(12) 



(13) 



E-/3 V(£;i-u)ni r 



is the partition function in the canonical ensemble. In the 
last expression the sum is taken over all possible states, 
|{ni}) of the quantum dot. To calculate Z{Ne) explic- 
itly we need to assume that energy levels of the dot are 
equidistant: 



(14) 



where A^e — Ni corresponds to the equilibrium number of 
the excess electrons (^ = /i); Z^xc is the partition func- 
tion of the thermal excitations. Now, let us substitute 
Eq. (O in Eq. (O: 



Zf_c 

_ 1 „-l3\U(N)+N^5E/2] 



(15) 
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where 



-f3[U(N)+N'^SE/2] 



(16) 



A'=~oo 



Here, we have extended one hmit of the sum to infinity 
since the Fermi energy is the largest energy scale of the 
problem. 

Thus, we are well-equipped to calculate the conduc- 
tance in the case of equidistant energy levels in the quan- 
tum dot at arbitrary ratio S. The widths of energy lev- 
els, r' and r[, in the quantum dot are energy depen- 
dent, random quantities. Let us assume that quantum 
dot is weakly coupled to the leads via multichannel tun- 
nel junctions: G'''^ = gY Nch <C e^//i, where Nch is 
the number of channels; GY is the conductance of one 
channel. Experimentally, this situation corresponds to 
the metallic grain coupled to the leads via oxide tunnel 
barriers iA This setup allows one to decrease fluctuations 
of the energy levels' widths, T\ and r[, by a factor of 
\/Nch- We also assume that the widths are slowly chang- 
ing functions of the energy, Ei. Then, F' and F[ can be 
considered constants and evaluated at the chemical po- 



tential: F: 



F^ rr 



— ^ , - 1 . - = F"". There is a simple 
relation between these widths and conductances of the 
corresponding junctions. In the case of spinless fcrmions: 



Let us substitute Eq. in Eq. Q: 



(17) 



G = ^, ^^2^< 



-(j[U(N)+N'^SE/2\ 



Y,Fe,{E,m 



N 



•K{l~nF[E,~ ^i + U{N)-U{N (18) 

To take advantage of the expression for occupation num- 
bers, Eq. (PJ, we need to map the sum over i onto the sum 
over j. As illustrated in Fig.|31 the mapping rule depends 
on the total number of excess electrons, N (compare with 
Eq. idni) written for iV = 0): 



5E + N6E, 



Feq{E,\Ne) = nj, 

E,-^,^E,-C + iC ^0 = (j-l 

where = C{^) is the energy of the highest occupied 
energy level in the dot with N excess electrons at T = 
plus 6E/2. We will also use the following identities: 

U{N) + = (^Ec + ^) - eVeN 

i?c + ^) (n-Ao + ^) +Ci, (19) 



where Ci = - (eVef /2{2Ec + SE); 

eVe 1 

Ao = + - 

2Ec + SE 2 



(20) 



has been chosen so that Aq = corresponds to the max- 
imum of the conductance peak; and 

U{N) - U{N - 1) = (27V - l)Ec ~ eVe 
^{2N-l)Ec-i2Ec + SE){Ao-l/2). (21) 

Substituting these results in Eq. (|18|l . we obtain 



G(Ao) ^ J_ 

Goo Dq 



— V e^^'o(w-Ao+i)' V- 

^0 V r 



where 



5-2(A'-Ao)eo + 1 ' 
(22) 



Do = exp{pCi)D'o = ^e-^«(^-^"+*)'; (23) 



N 



eo = P{Ec + SE/2); Goo = G^G'-fiG^ + G^ is the clas- 
sical, Ec,SE <^ T <Si n, limit of the conductance. We 
have also used the identity: 1 — npiE) = (1 -I- e~^^)~^. 
Formula (|22(l is the general expression for the linear con- 
ductance in the spinless case for equidistant energy levels 
in the quantum dot at arbitrary values of Eq, SE and T. 

One can immediately prove the following properties of 
the conductance First of all, G(Ao) = G{Ao+M), 
where M is an integer. In the gate voltage units, AVe — 
{2Ec + SE)/e is a period of the conductance oscillations. 
This property reflects symmetry with respect to adding 
(removing) an electron to the quantum dot. Secondly, 
due to the electron-hole symmetry, conductance is an 
even function of Aq: G(Ao) — G(— Aq). 

The linear conductance H22|) as a function of the di- 
mensionless gate voltage, Aq, at SE = Ec is plotted in 
Fig. 21 for different temperatures. At low temperatures 
there are sharp Coulomb blockade peaks. At high tem- 
peratures, T Ec,SE, Coulomb blockade is hfted and 
small oscillations of the conductance can be observed. 
These oscillations are slightly non-sinusoidal and given 



O 




FIG. 4: Coulomb blockade oscillations of the conductance 
as a function of the dimensionless gate voltage, Aq at 6E = 
Ec- Curves are plotted for different temperatures: T/Ec = 
T/SE = 1/5 = 0.2, 0.35, 0.5, 0.7, 1, 2, 5. 
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by the following asymptotic formula: 



G(Ao) - G(Ao) 



-£o/4 r 



cos(27rAo) 



Goo -v/Eo 

+ e~2^'/"° cos(47rAo) + 0(e-3^'/"°)l (24) 



where G(Ao) is the average value of the conductance. 
The second term in Eq. (|24l) is due to inherently non- 
sinusoidal nature of the conductance oscillations, see 
Fig. ^ To derive this expression one can use Poisson's 
summation formula. 

To study the line shape of a separate peak let us con- 
sider the limit of large charging energy: Ec ^ T, SE 
or, equivalently, Eq S> 1,(5 in the dimensionless units. In 
Eq. (|22|l only N — —1,0 terms in the sum over N give 
substantial contribution to the conductance near Aq — 0; 
all other terms are exponentially suppressed. Besides, 
sum over j at the iV = — 1 is 0(e"^^") and, therefore, 
can also be neglected. Hence, line shape of the conduc- 
tance peak at Aq = is given by 



G(Ao) 

Goo 



1 + e-2eoAo 1 



,-j<5+2eoAo ■ 



(25) 



It is more instructive to rewrite this equation as follows 
G(0o) 5 ^ n 



Go 



j 



(26) 



where (j)o = — vj^^] v}^^ is chosen so that c/jq — 
corresponds to center of the conductance peak. In the 
classical regime, T ^ SE, line shape of the conductance 
peak is given bji 



G(0o) 

Goo 



2 sinh(/3e0o) 



(27) 




FIG. 5: The exact line shape of the conductance peak, 
Eq. plotted at 1/5 = T/SE = 0.4 (solid curve). Dotted 

curve corresponds to blindly applying classical regime formula 
(I27I I at T/SE = 0.4. Dashed curve corresponds to blindly ap- 
plying formula (I28II at T/SE — 0.4. The argument of the plot 
is linear function of the gate voltage, (po- 



In the opposite limit of SE ^ T: 



Goo 



2[1 + cosh(/3e(/)o)]' 



(28) 



The exact line shape of the conductance peak, Eq. H26|l . 
at T/SE — 0.4 is shown in Fig. \^ On the same figure 
we also plotted two conductance peaks in the limiting 
cases, Eqs. (|27|l and (|28|l . out of their validity region at 
T/SE — 0.4. Nevertheless, it is interesting that the ex- 
act conductance peak is higher than both of the limiting 
cases peaks. The peak's height is given by 



G(0) 

Goo 



2 ^ 1 

3 



1/2, 
(5/4, 



T-» SE 

SE:$>T:$>r, 



(29) 

Temperature dependence of the conductance peak's 
height was numerically calculated in the Ref. UXl see 
Fig. 2 there. 



IV. APPLICATION TO TUNNELING 
THROUGH QUANTUM HALL EDGE STATES IN 
A QUANTUM DOT 

Formulas for the linear conductance in the case of 
equidistant energy levels in a dot and spinless fermions 
derived in the previous section can be applied to a num- 
ber of physical problems. 

Let us consider, for example, a quantum dot formed by 
confining a two-dimensional electron gas by a circularly 
symmetric electrostatic potential, U{r). We assume that 
U{r) is zero at the origin and takes large value ai r — R, 
where R is the radius of the dot, Fig. Let us apply 
a strong magnetic field, B, perpendicular to the plane of 
the dot. This situation corresponds to the quantum Hall 
regime and was reviewed in Ref. 

To solve the one-electron Schrodinger equation in this 
geometry it is convenient to choose the symmetric gauge. 
Then, angular momentum is, clearly, an integral of mo- 
tion. In each Landau level, n, states with larger angu- 





a) 



b) 



FIG. 6: Geometry and spectrum of the quantum dot states: 
a) symmetric gauge eigenstates, |m), including edge state are 
shown schematically; b) energy spectrum of the eigenstates 
with different angular momenta, two lowest Landau levels are 
shown. 
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lar momentum, |m), are localized further from the ori- 
gin, near a circle with radius Rm = Ih \/2(to + 1), where 
Ih = \fhcfeB is the magnetic length, c is the speed of 
light. The presence of the confinement potential leads 
to an increase in energy for the symmetric gauge eigen- 
states with i?m of order or larger than i? (or, m of order 
or larger than mp-, see Fig. |H|d). For the states with 
Rra ~ i? an electron is influenced by both the electric 
field of the boundary, E{K) — U'{R)/e, and strong, per- 
pendicular to the electric, magnetic field. Thus, near the 
edge electron executes rapid cyclotron orbits centered on 
a point that slowly drifts in the direction of E x B, that 
is, along the boundary. Thus, Quantum Hall edge states 
are formed, Fig.^Ji. It is important to notice that in this 
closed geometry electron system has only one edge. In 
this consideration we also assume that Ih ^ ^• 

For simplicity let us consider the case when only the 
zeroth Landau level crosses the chemical potential, that 
is, there is only one type of edge states. This corresponds 
to a sufficiently strong magnetic field so that filling factor, 
V is equal to 1, Fig. EId. 

Now, we are ready to consider transport through this 
type of quantum dot in the strong magnetic field. Let us 
weakly couple it to two leads and apply an infinitesimally 
small bias voltage between them. An electron from the 
left lead can now tunnel into dot's edge state and then 
tunnel into the right lead as illustrated by dashed lines 
in Fig.EJ^. 

The energy spectrum of edge states can be linearized 
as follows 



Cm = /i + 

am 



(m — mp) . 



(30) 



Thus, energy levels of the edge states are equally spaced 
with the spacingi^ 



5E 



dm 



dRn 



dRr, 



dm 



eE{R)^. (31) 



This fact makes formulas derived in the previous section 
applicable to this problem. Essential assumption here 
is that dispersion curve. Fig. is almost linear in the 
range of angular momentums: \m — mp\ ^ Am, where 
Am = ma.x{l,T/SE). 

In the case at hand, spacing, SE is inversely propor- 
tional to the size of a dot just like charging energy, Ec- 
Hence, their ratio does not depend on the size of a dot 
and is given by 



5E 

E^ 



e E{R) 
a B '■ 



(32) 



where e is the dielectric constant of the media around 
the interface; a = /Tic is the fine structure constant. 
Therefore, in this case oscillations of the conductance 
given by Eq. H22|l are determined by only one parameter 
5 = 5E/T. 

In conclusion, let us consider the case of an arbitrary 
shaped quantum dot. In this case, m is just the index of 



an edge state and no longer associated with the angular 
momentum. The phase along the boundary for the m-th 
edge state is 



dx fcjTi(x), 



(33) 



where km{x) is the corresponding wave vector, x 
parametrizes the boundary, and L is its length. The 
phase difference between two consecutive edge states is 



27r 



dx [fcm+i [x) - km (x)] , (34) 



where km+i{x) - km{x) = {£^+1 - £m) /Tiv{x), v{x) = 
cE{x)/B is a drift speed along the boundary. Then, the 
spacing between edge states' energy levels is 



5E = 2TTh 



dx 


-1 


dx 


-I 


Jo ^i^) 




Jo E{x)_ 





(35) 



Though the electric field E{x) at the boundary slightly 
changes as one goes from one edge state to the other, this 
effect is small and we neglect it. Therefore, the energy 
levels of the edge states are equidistant with the spacing 
given by Eq. JHSJ. 

In the case of the circularly symmetric quantum dot, 
E{x) is constant, and one can easily perform the integra- 
tion in Eq. (^5)) . This leads to the previously obtained 
expression for the level spacing, Eq. H31(l . 



V. LINEAR CONDUCTANCE IN THE SPIN-i 
CASE 

Formula Q for the linear conductance in the spinless 
case can be easily generalized to the spin-i case by count- 
ing each energy level twicet^ii 

G = 2—5]-^ ^ P^,{N^)F^,{E,,\N^) 

^[l-npiEi- ^ + U{N)-U{N -I))], (36) 

where Feq{Ei^ |iVe) is the occupation number of the quan- 
tum dot's energy level i with a spin- up electron, («, t) in 
the canonical ensemble: number of electrons in the dot, 
Ne, is fixed. 

As in the spinless case, to carry out analytical con- 
sideration we have to assume that energy levels in the 
quantum dot are equally spaced. Presence of the spin 
degeneracy makes the calculations more complicated. 

First of all, let us find the occupation number 
Feq{Eii \Ne). Lct US cousidcr spin-up and spin-down elec- 
tron subsystems. Groimd state energy of the system is 



E„ 



^1 



6E 



5E 
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where = {Ng)a — Ni/2 is the number of excess 
electrons in the spin-cr subsystem; Ni is chosen even. 
N^^ ^ Ni = N \s the total number of excess electrons 
in the quantum dot; [N^ — Ni)/2 = Sz is z-component 
of the total electron spin. Using these identities, one can 
find that 



1 



E„ = -N^SE + Sl5E. 



(37) 



While Sz is subjected to the thermodynamic fluctuations, 
N is fixed. 

The occupation number in question, _Feg(i?i| |A^e) = 
riji, is known if, in addition to N, the z-component of 
the total spin, 5*2, is fixed. In the case of an even num- 
ber of electrons, parameter C for the spin-up electron 
subsystem is equal to C, + Sz5E, given Sz^ see Figs. [7Ji 
and[7|3. Occupation numbers in the spin-up subsystem 
at fixed Sz are given by Eq. © with the appropriately 
chosen parameter C,': nj-s^- 

However, z-component of the total spin, Sz, is not fixed 
but subjected to the thermodynamic fluctuations. There- 
fore, to find the occupation numbers, n||, we have to 
account for all possible values of Sz- 



where 



-Sis 



NJ2 



(39) 



is the probability that z-component of the total spin of 
the quantum dot is equal to Sz- Substituting Eqs. H39(l 
and © in Eq. we obtain 



oo ^ 



We/2 



Since number of electrons in the quantum dot is even, Sz 
may take only integer values. Therefore, 

oo „ 

°° r 2 ^ ^ 
nf^ = ^(-l)'-ie-L^+(^'-*H* , 



a) 



Nc/2 

J2 P{Sz)n,.s., 
= -Nj2 



j=2 
j=l 



C + 5E 



b) 



(38) 



j=2 
j=l 
j=0 
j=-l 



where we extended limits of the sum over Sz to infinities 
since the Fermi energy is the largest energy scale in the 
problem. Separating m = 2r — 1 and m = 2r parts 
of the sum, where r is a positive integer, we obtain final 
expression for the occupation numbers in the case of even 
number of electrons: 



where 



r—l r—1 



-r{r+2j-l)5 



A{6) = 



— — OO 



(40) 



(41) 



c) 



j=2 
j=l 
j=0 
j=-l 



? + (S -1/2)5 E 



d) 



j=2 
j=l 
j=0 
j=-l 



FIG. 7: Parameter C,' of the spin-up electron subsystem for 
even number of electrons: a) at = 0, b) at arbitrary integer 
Sz {Sz = 1 case is shown); for odd number of electrons: c) 
at Sz ~ 1/2, d) at arbitrary half-integer Sz {Sz = 3/2 case is 
shown) . 



In two limiting cases 

26-''/'' [1 + (e-*')] , 5E^T 



A{5) = 



T-»5E 



Analytical expression for the high-temperature limit of 
A{S) can be obtained using Poisson's summation formula. 

One can easily prove the following properties of the 
occupation numbers valid at arbitrary temperature: 



= 1-n 



A{S)e^/\ 



(42) 
(43) 



They are valid due to the electron-hole symmetry and 
similar to the following properties of the Fermi-Dirac 
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distribution: npiE) = 1 — np{—E) and e^^npiE) 



e'^^npi-E) = 1. 

Similarly, one can find occupation numbers in the case 
of odd number of electrons in the quantum dot, N^- En- 
ergy level which contains one electron at T = will be 
referred to as j = level. In this case electron-hole sym- 
metry corresponds to j — s- — j transformation. Parameter 
C of the spin-up electron subsystem at a given 5*2 is equal 
to C + (5*2 - 1/2)5^;, see Figs.|7|: and|7|l. Therefore, 

Afe/2 



We/2 , x2 



We/2 



.(44) 



Since number of electrons in the quantum dot is odd, 
may take only half-integer values. Separating odd and 
even parts of the sum over m, we obtain: 



^ oo oo 

A{5) ^ ^ 



-r{r+2j)S 



(45) 



Property of the electron-hole symmetry reads as follows: 

"?f = l-<T- 

It turns out that there exists simple relation between 

and n°'^ occupation numbers: 



(46) 



This property is the analog of np{E) = e~^^np{—E) 
one of the Fermi-Dirac distribution. It will allow us to 
get rid of n°'^ occupation numbers in the final expression 
for the conductance. 

Now we are in a position to find the probability that 
a dot, in thermodynamic equilibrium with the reservoirs, 
contains electrons, Peq{Ne). Equation H13|l written for 
the spinless case is still applicable if we keep in mind that 
energy levels in the quantum dot are doubly degenerate. 
Partition function of the dot's internal degrees of freedom 
in the canonical ensemble is 

(47) 



where 



{".} 



(48) 



is the partition function of the spin-up electron subsys- 
tem in the canonical ensemble. Mathematically, expres- 
sion for Z{NJ) is identical to the one for Z{Ne) in the 



spinless case. Thus, one can directly apply the result 
obtained previously, Eq. (|14|l : 



ZiNj) 



(49) 



where iVj = Nijl is the equilibrium number of electrons 
in the spin-up subsystem. Similar result is valid for the 
partition function of the spin-down electron subsystem, 
Z(N^). Substituting these results in Eq. (|47|l we obtain: 



Z(7Ve) = Z, 



2 



We We 



We 



Afj=0 ATe -O 

The exponent can be simplified as follows 



hence, 



Z{N,) 

y2 

exc 



iVe We 



E E e"(^= ''/'*(5 



Nl+Nt,N^ 



We We/2 

^ g-W"5/4 ^ g-(WeT-We/2)^5 ^ g-W25/4 ^sH 



wJ=o 



E 

S = -We/2 



where in the second equality we took advantage of the 
delta symbol. Sum in the last line is taken over integer 
values of s if is even or half-integer values of s if is 
odd. Limits of the sum over s can be extended to infini- 
ties since we assume that ^ ^ T. According to Eq. 
probability that quantum dot, in thermodynamic equi- 
librium with the reservoirs, contains Np electrons is 



s=-A'e/2 



We/2 



We/2 

-f3U(N) g-s^5 

S = -We/2 



N=ev s N=od s 

where U{N) = U{N) + N'^5E/4; and we used the fact 
that TV and have the same parity since Ni is chosen 
even. Therefore, sums over N — ev and N — od are 
taken over = 0, ±2, ±4, . . . and N = ±1, ±3, . . . values, 
respectively. At this point in the calculation we need to 
specify whether the total number of electrons in the dot 
is even or odd: 

P (TV)-/^^')"'^""''^''^'- (50) 
i^e,(iVej - < A{5)e-P^^^\ N is odd ' ^^"^ 



where 



E 



-f3U{N) 



MS) E 



-/3C/(W) 



(51) 



N=odd 
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Now we are prepared to calculate the conductance, 
Eq. in the case of the equidistant double degen- 

erate energy levels in the dot at an arbitrary 5E/T and 
Ec/SE ratios. Similarly to the consideration in the spin- 
less case we assume that quantum dot is weakly coupled 
to the leads via multichannel tunnel junctions, and tun- 
neling widths of the energy levels in the quantum dot, T\ 
and r[, are slowly changing functions of the energy, Ei. 
Then, these tunneling widths can be considered constants 
and evaluated at the chemical potential: « = F'; 
rr ~ = F*". Furthermore, they can be expressed via 
conductances of the corresponding junctions: 



hG' 5E 



hC SE 



(52) 



There is an additional factor of 1/2 here compared to 
the spinless case, Eq. ((TTIl . due to the double degeneracy 



of each energy level in the quantum dot. First of all, 
let us break the sum over Nf, in Eq. (|36|l in two parts: 
N — even and N — odd, and apply Eqs. and (|^ : 



. j ^ r r 

t I j = Q j = 2 j = i 

t I j = p j = i j = o 

t \ T'T" j = o = 

A : A : A 

V V V 

a) b) c) 

FIG. 8: Mapping of the sum over i onto the sum over j for 
a) A'' = 0; b) even A'' (A'^ — —4 case is shown); c) odd A*' 
{N — —3 case is shown). 



G 



G^G'' S 



E 



-I3U{N) 



J2 Feq{E^^ \Ne) {1 -nF[E,-iJL^ U{N) - U{N - 1)]} 



\N—eve7i 



+A{5) J2 e-'^^W^Fe,(i?,T|iVe){l-n;^[i?,-M + f/W-f/(A^-l)]} 

N=odd i=l J 



(53) 



r 



To take advantage of the occupation numbers we derived, 
Eqs. I|4U|) and H45|) . we need to map each of the sums over 
i onto the sum over j. The first sum over i in Eq. I|53|) 
is taken at even number of excess electrons, see Figs. 
andlSb, hence 



E,-^l = E,-C + {C -0 = (j- 



i j 6E+^NdE. 



Remember that by properly choosing "zero" of the gate 
voltage we put C — ^. The second sum over i is taken 
at odd number of excess electrons, see Figs. and 
therefore 



where C2 = -ieVe)'^ j [\Ec + 5E); 



A = 



eVe 



C2, 



(54) 



2{Ec + 8EI\) 2 
is the dimensionless gate voltage, A = corresponds to 
a position of the conductance peak; and 



E,-ix = E,~C: ^{C -0 = jSE + ^NSE. 
We will also use the following identities: 



U{N) - U{N - 1) = (2A^ - l)Ec - eVe 
= (2A^ - l)Ec -2(ec + -^^ - ^ 



U{N) = [Ec + ^]N'^- eKN 



Substituting these results in Eq. (|53|1 we obtain 



G(A) 

Goo 



S 
D 



N—even 



-e(A^-A+l/2)^ 



E ' 



1 



-U-l/i)S-2{N-A)e I 



-0-l/4)5+2(Ar+l-A)e _^ i 



(55) 
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where 

D^expipC2)D' ^ J2 e-'^^-^+'/^^" + A{S) ^ ^-e{N-A+i/2f . (gg) 

N—even N—odd 

I 



e^pl^Ec+^-^y, (57) 

Goo = G^C/iG' + G"^) is the high-temperature, 
Ec, SE ^ T ^ /i, Umit of the conductance. To ehminate 
^od fpQjjj ^jjg final expression we used useful property of 
the occupation numbers given by Eq. 14t)|l . 

Formula (|55|l is the main result of this paper. It is the 
analytical expression for the linear conductance in the 
spin-i case for equidistant energy levels in the quantum 
dot. One can use formula (|55|) to plot Coulomb block- 
ade oscillations of the conductance as a function of the 
dimensionless gate voltage, A, at arbitrary values of Ec, 
6E and T. Particularly, when all energy scales are of the 
same order: Ec ^ SE ^ T, numerical calculation is a 
breeze, Fig. 

One can immediately notice the following properties 
of the linear conductance. First of all, G(A) = G(A + 
2M), where M is an integer. In other words, {4Ec + 
SE) /e is the conductance period in the gate voltage units. 
This property reflects symmetry with respect to adding 
(removing) two electrons to the quantum dot. Secondly, 
conductance is a symmetric function with respect to the 
center of a valley. A' = M + 1/2, where M is an integer. 
That is, G(A - A') = G(A' - A). This is a reflection 
of the electron-hole symmetry. These properties of the 
conductance oscillations are not generic. They are valid 
due to the assumption of equally spaced energy levels in 
a quantum dot. 

At high temperatures, T ^ Ec,SE, conductance 
peaks overlap and their maximums become almost 



1- 




FIG. 9: Linear conductance oscillations as a function of the 
dimensionless gate voltage. A, at SE — Ec- Curves are 
plotted for different temperatures: T/Ec ~ T/5E = 1/5 = 
0.15, 0.3, 0.5, 0.7, 1, 2, 5 using Eq. iCT . 



equidistant, Fig.|^ As a result, instead of separate peaks, 
the conductance in this limit has oscillatory behavior. 

Let us find temperature dependence of the conduc- 
tance in the valleys. In the sequential tunneling approx- 
imation the conductance in the valleys decays exponen- 
tially as T ^ 0, Fig.^l At low temperature number of 
electrons in a dot in the valleys is almost quantized. We 
will call the valley "odd" ("even") if it corresponds to 
odd (even) number of electrons in the dot. We find that 
at any temperature the conductance in the odd valley is 
larger than that in the even one. Fig. ^] This feature is 
robust with respect to the distribution of energy levels in 
a quantum dot. 

However, it is important to mention that at low tem- 
peratures, T < Tin, where 




cotunneling^^ will dominate sequential tunneling contri- 
bution to the conductance in the valleys. Therefore, tem- 
perature dependence of the conductance in the valleys. 
Fig. EH is valid only for the temperatures T > Tin. 

Let us analyze the limit of large charging energy: 
Ec 3> T,SE or, equivalently, e ^ 1,(5 in the dimen- 
sionless units. In this limit, two adjacent peaks in the 
conductance have exponentially small, ~ e"^'^/"^, over- 
lap with each other. Thus, it makes perfect sense to 
study the line shape of a separate peak. Let us deter- 
mine line shape of the conductance peak near A = 0. In 
the numerator of Eq. (|55() only the = term in the 
sum over A^ survives; moreover, at A^ = second term in 




0.2 0.4 0.6 0.8 



T/8E 

FIG. 10: Temperature dependence of the conductance in the 
odd, G(-l/2)/Goc. and even, G(l/2)/Goo valleys at 5E = Ec- 
Upper curve corresponds to the conductance in the odd valley. 
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square brackets is 0(e~^^). In the denominator, only the 
N — —1,0 terms matter. Hence, the line shape of the 
conductance peak near A = at arbitrary SE/T ratio is 
given by 



1 



^(,5)e-/3e0 ^ 1 + e-(j-l/4)5+/3e.^ 
j = -oo 



(59) 



where we used the following identity: 

2eA = (3e (Ve ~ = (3e<j>. 

Clearly, cf) = corresponds to A = 0. In the classical 
regime, T ^ SE, the line shape of the conductance peak 
is given byi 



Goo 



2 sinh(/3e0) ' 



(60) 



Formally, this equation is identical to that of the spinless 
case, Eq. (|27f) . Nonetheless, the values of Goo are differ- 
ent in these two cases by a factor of 2. This is due to 
spin degeneracy of each energy level in the spin-i case, 
compare Eqs. and lO- In the hmit of 6E > T& 



GW 



Goo 3 + 2V2cosh(/3e(/)+ i^- iln2)' 
This peak has its maximum at 

e(l)LT = e(b,n (T « 6E) = -\sE + i^T, 



(61) 



(62) 



and is symmetric with respect to this value: G(0 — 
4>lt) — G{(I>LT — 0)- The exact line shape of the con- 
ductance peak, Eq. (|59|l . at T/SE — 0.4 and two limit- 
ing cases conductance peaks, Eqs. and (jHU, plotted 




FIG. 11: The exact line shape of the conductance peak, 
Eq. lEHJ, plotted at 1/5 = T/SE = 0.4 (solid curve). Dot- 
ted curve corresponds to blindly applying classical regime for- 
mula, Eq. at T/6E = 0.4. Dashed curve corresponds to 
blindly applying formula 1611 at T/SE = 0.4. The argument 
of the plot is linear function of the gate voltage, <j)- 



out of their validity region at T/5E = 0.4 are shown in 
Fig. ^2 As in the spinless case, the exact conductance 
peak is higher than both of the limiting cases peaks. 

The position of the peak's maximum, (jjm — 4>m{T)^ 
is shifted to the left from its high temperature limit, 
(l>CL = 4>m{T > SE) = 0. It is determined by the 
equation for (/>,„: G' (</).,„) = 0, where G(0) is given by 
Eq. H59|l . The dimensionless position of the peak's max- 
imum, e(j)rn{T)/5E, as a function of the temperature, 
T/5E, is numerically plotted in Fig. [T^ 

The conductance peak height is Gmax — G[4>m)- In 
the limiting cases: 



Gr, 



Go 



G{4>ra) 
Goo 



1/2, T->5E 
(3-2^/2)(5, (5£;>r>r, 



(63) 

Peak's height as a function of the temperature, T/5E^ 
can be plotted numerically. Fig. ^] 



VI. CONCLUSIONS 

We have studied Coulomb blockade oscillations of the 
linear conductance through a quantum dot weakly cou- 
pled to the leads via multichannel tunnel junctions in the 
sequential tunneling approximation. To obtain analyti- 
cal results we have assumed that the energy levels in the 
dot are equally spaced. The electron-electron interaction 
in a quantum dot has been described by the constant 
interaction model; though, thermal excitations with all 
possible spins have been taken into account. 

The linear conductance in the spinless case is given by 
Eq. H22|) . It is valid at arbitrary values of Ec, SE and 
T. The line shape of an individual conductance peak at 
arbitrary ratio S = SE/T is given by Eq. H26|l . Exact con- 
ductance peak is higher than both of the limiting cases 
peaks at any gate voltage as is illustrated in Fig|5| An 
analytical expression for the height of the conductance 
peak at any ratio S is obtained, Eq. (|29|l . 



-0.05 - 



0-1 
CO 



-e- 



-0.1 




-0.15 



-0.25, 



FIG. 12: Temperature dependence of the dimensionless 
peak's maximum position, e4>m{T) /SE. In the low tempera- 
ture limit: e4>m{0)/5E = — 1/4 according to Eq. 
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In Sec. IV we applied the spinless case theory result 
to the problem of the transport via a dot in the quan- 
tum Hall regime. Energy levels in a dot in this case are 
equidistant with the spacing given by Eq. (|35ll . 

Linear conductance in the case of spin-^ electrons at 
arbitrary values of Ec, 6E, and T is given by Eq. (|55|l . 
In particular, this equation allows one to plot the con- 
ductance oscillations in the regime when the charging 
energy, level spacing in the dot, and the temperature are 
all of the same order. Fig. |^ We find that the period 
of Coulomb blockade oscillations is doubled compared to 
the model with a continuous electronic spectrum in the 
dot. Equation is the main result of the paper. 

We also find that conductance in the odd valley is 
larger than that in the even one at any temperature. 
Fig. ^1 The difference between conductances has the 
largest value at T w 0.3 6E (at SE = Ec). The sign 
of the difference is the same as for the quantum dot in 
the Kondo regimer^' Kondo effect takes place at very low 
temperatures, T < Tk <C Tin, where Tk is the Kondo 
temperature, and leads to the logarithmic enhancement 
of the conductance in the odd valley. Our consideration 
shows that even-odd asymmetry exists at much higher 
temperatures. 

Line shape of the conductance peak is given by 
Eq. H59|l . As in the spinless case, the conductance peak is 
higher than both of the limiting cases peaks at any gate 
voltage, Fig^] As we increase the temperature peaks' 
maximums shift and become more equidistant. Fig 1121 
The peak's height as a function of the temperature is 
calculated numerically and plotted in Fig. E| 

Though we have found physical system which has 
equidistant energy levels in the spinless case, see Section 
IV, we are not aware of any such system in the spin-^ 
case. In the case of a chaotic quantum dot Wigner-Dyson 
model gives a fairly good approximation for the distribu- 
tion of the energy levels of the dot. If we had assumed 
Wigner-Dyson distribution of the quantum dot's energy 




0.5- 



0.2 0.4 0.6 

T/5E 

FIG. 13: Height of the conductance peak, Gmax/Goo, as 
a function of the temperature, T/SE. Dashed curves corre- 
spond to two limiting cases peak heights, Eq. 15311 . 



levels then we would have had to give up the hope of 
finding a solution. It goes back to the very difficult prob- 
lem of finding occupation numbers of the dot's energy 
levels in the canonical ensemble. The only way to solve 
it is to assume that energy levels in the quantum dot are 
equally spaced. Then one can use the bosonization tech- 
nique to find the occupation numbers. Assumption of the 
equidistant energy levels is in line with the level repulsion 
property of the Wigner-Dyson distribution. Therefore, 
the analytical consideration of this reasonably simplified 
model, in our opinion, is a significant step forward in the 
solution of the general problem. 

Though we do not expect our quantitative results to 
precisely describe a quantum dot with random energy 
levels, they certainly give correct order of magnitude for 
the conductance oscillations and their generic features. 
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APPENDIX: OCCUPATION NUMBERS IN THE 
CANONICAL ENSEMBLE 

In the canonical ensemble, where the total number of 
particles, TVg, is fixed, it is difficult to calculate fermionic 
occupation numbers, Eq. (jS]), directly. This is true even 
in the case of equidistant energy levels separation. 

Fortunately, the bosonization technique allows one to 
express fermionic field annihilation and creation oper- 
ators, "0 ^-iid "0^1 in terms of the bosonic annihilation 
and creation, a^'s and aj's, and ladder, U and J7^, 
operatorsii^ 

^ ^e~''^''^^e-*^^(^)C/te-^x(^), (A.l) 

where 

x'{-)^"-^N + ^Y,J^e-^^^al■ (A.2) 

L is the length of the artificial system; fcj? is the Fermi 
wave vector; N = Ne — Ni is the number of excess elec- 
trons operator; q is the wave vector; and x is the coor- 
dinate. Since these bosons naturally exist in the grand 
canonical ensemble, one can fix total number of excess 
electrons, TV, and calculate occupation numbers in the 
bosonic basis, Eq. 0. 
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Occupation numbers, Eq. have the following prop- 
erties: 



in the high temperature limit, T ^ SE: 



(A.3) 



They are similar to np{E) — 1 — np{—E) and 
nF{E)el^^/'^ = nFi-E)e'^^/^ ones of the Fermi-Dirac 
distribution and reflect the electron-hole symmetry. 
Combining these two properties one can get the recur- 
sion relation: 



rij+i 



1 



(A.4) 



In the limit of low temperature, SE ^ T, for j > we 
obtain 



(3i+3)5 



(A.5) 



5 sinh [(j 



' cosh'^ [ ( j 



0{S')iA.6) 



thus, we find first correction to the Fermi-Dirac distri- 
bution. The form of the correction remains valid even 
for the slightly non-equidistant energy levels in the dot. 
If this is the case we need to define ratio S in the last 
expression as 5 — SE/T, where SE is the mean level 
spacing. 
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